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On the basis of a coherent state representation of quantum noise operator and an ensemble 
averaging procedure a scheme for quantum Brownian motion has been proposed recently [Banerjee 
et al, Phys. Rev. E 65, 021109 (2002); 66, 051105 (2002)]. We extend this approach to formulate 
reactive flux theory in terms of quantum phase space distribution functions and to derive a time 
dependent quantum transmission coefficient - a quantum analogue of classical Kramers'-Grote-Hynes 
coefficient in the spirit of Kohen and Tannor's classical formulation. The theory is valid for arbitrary 
noise correlation and temperature. The specific forms of this coefficient in the Markovian as well as in 
the non-Markovian limits have been worked out in detail for intermediate to strong damping regime 
with an analysis of quantum effects. While the classical transmission coefficient is independent of 
temperature, its quantum counterpart has significant temperature dependence particularly in the 
low temperature regime. 

PACS numbers: 05.40.-a, 82.20.-w, 82.20.Db, 02.50.Ey 



I. INTRODUCTION 



The classic work of Kramers 1 on the diffusion model of chemical reactions in 1940 forms the dynamical basis of 
modern rate theory of activated processes. With the development of ultrafast lasers and time resolved detection 
techniques since mid seventies it has been possible to monitor the details of the pathways of a reaction on a micro- 
scopic scaled. This has provided a major impetus and also generated a renewal of interest in theoretical development 
in chemical dynamics, particularly in the problem of reactions in condensed media. Over the decades the field has 
grown in various new directions, e.g. , extension of Kramers results to non-Markovian regim e 4 i 5 i 6 , generalizations to 
higher dimensions 7,8 , inclusion of complex potentials' -, generalization to open systems 1 "^, analysis of semiclas- 
sical and quantum effecta 14 i 15 i 16 i 17 i 18 i 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 i 29 , thermal ratchet and molecular motors^ etc. These 
developments have been the subject of several reviews and monographs. We refer to 1 " 7 ^. 

It is now well known that the rate constant of a chemical reaction may be realized as a typical transport coefficient 
in phase space. This transport coefficient can be expressed in the form of a Green-Kubo formulae for the flux-flux 
auto correlation function calculated at the transition gtateiLSSiSiSii&SiiSI which acts as a dividing surface between 
the reactant and the product states. The theory of reactive flux has been developed both in classical and quantum 
mechanical contexts to obtain the rate coefficient of a chemical reaction. Recently Kohen and Tannor— have developed 
a method to calculate the classical rate coefficient and hence the transmission coefficient as a function of time in 
terms of phase space distribution function in the high friction regime. The method is simple and directly addresses 
the relevant noncquilibrium dynamics at the barrier top. The theor y 36 : 37 not only recovers the Kramers-Grote-Hynes 
transmission coefficient in the long time or asymptotic limit but also traces the nature of its transient dynamics and 
its approach towards the steady state. 

The key element of the Kohen- Tannor formulation— is that it allows us to express the correlation function with the 
help of a probability distribution function which describes dynamics of the reaction co-ordinate undergoing Brownian 
motion in phase space at the barrier top without resorting to a solution of the generalized Fokker-Planck equation. 
It is worthwhile in this context to ask whether the theory can be extended to quantum domain to include the generic 
quantum mechanical effects. The essential requirement for such a theory is obviously the quantum phase space 
distribution function. Very recently, based on a coherent state representation of quantum noise operator and an 
ensemble averaging procedure, we have developed 24 ' 25-26,27 a c-numbcr quantum generalized Langevin equation. The 
theory is valid for arbitrary temperature and noise correlation. We extend the scheme in the present context to make 
use of a true quantum probability distribution function of the relevant dynamical variables in c-numbers and express 
the rate constant or transmission coefficient accordingly. Our object is thus two-fold, i.e. , 

(i) to express the quantum mechanical rate constant in terms of the c-number quantum phase space functions. 
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(ii) to derive a time dependent quantum Kramers-Grote-Hynes (KGH) transmission coefficient which is valid for 
arbitrary noise correlation and temperature for intermediate to strong friction regime. 

Two pertinent points are now in order. First, generally speaking the phase space distribution functions used in 
quantum mechanical contexts, e.g. , in quantum optics are quasi-probability distribution functions like Wigner, 
or Glauber-Sudarshan functions 38 . They often become singular or negative in the fully quantum domain and can 
therefore not be regarded as true probability distribution functions as their classical counterparts. We stress that 
the present formalism is free from such difficulties. Second, the expression for the rate constant as a time integral 
over flux-flux auto correlation function has been evaluated earlier both analytically and numerically using path 
integral approaches. Although for equilibrium properties the imaginary time propagator has proved to be very 
much successful 23 , particularly for developing quantum transition state theory^, it is not easy to extract dynamical 
information for noncquilibrium problems^ 2 , because of the oscillatory nature of the real time propagator in many 
situations. The c-number phase space method that we use here, on the other hand, being independent of path 
integral approach is an alternative for calculation of the rate without these problems. 

The outlay of the paper is as follows; After a brief review of our recent scheme for quantum Brownian motion we 
introduce the conditional probability function of c-numbers at the barrier top and an equilibrium distribution function 
at the reactant well in Sec. II. In the following section the reactive flux expression for rate constant is formulated in terms 
of these c-number phase space functions. The time dependent general quantum Kramers-Grote-Hynes transmission 
coefficient has been derived in the spirit of Kohcn and Tannor formalism. We work out in detail the Markovian and 
non-Markovian limits of this expression. The paper is concluded in Sec. IV. 

II. QUANTUM BROWNIAN MOTION ; EVOLUTION OF C-NUMBER PHASE SPACE FUNCTIONS 



The kinetics of a particle undergoing a chemical reaction in a medium can be described by the following generic 
system-reservoir Hamiltonian 



Here q and p are the co-ordinate (here considered as reaction co-ordinate) and momentum operators of the particle, 
respectively and the set {qj,Pj} is the set of co-ordinate and momentum operators for the reservoir oscillators linearly 
coupled to the system through coupling constant Kj. V(q) denotes the external potential field which, in general, is 
nonlinear. The Hamiltonian contains counter-terms quadratic in q to ensure the independence of barrier height and 
coupling constants. The operators follow the usual commutation relations [q,p] — ih and [<2j,Pj] = iHSij. Over the last 
several decades the model and its variants have widely served as the fundamental paradigm for quantum Brownian 
motion in chemical dynamics, condensed matter physics, quantum optics and allied areas within various approximate 
descriptions. Our primary aim in this section is to search for a true c-number phase space function which describes 
the evolution of quantum Brownian motion of the particle at the barrier top and its thermalization in the reactant 
well. 

Eliminating the reservoir degrees of freedom in the usual way we obtain the operator Langevin equation for the 
particle 



A. Quantum Langevin equation in c-numbers 




(1) 




(2) 



where the noise operator F(t) and the memory kernel j(t) are given by 




(3) 



j 



and 




Kj COSOJjt 



j 
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or in the continuum limit 

/>oo 

7(f) = / k(lu)p(lu) cos cot duj (4) 
Jo 

with Kj — io 2 y The masses have been assumed to be unity. p(u>) represents the density of the reservoir modes. 

Eq.J2J) is an exact Langevin equation for which the noise properties of F(t) can be defined using a suitable canonical 
initial distribution of bath co-ordinates and momentum. Our first task is to replace it by an equivalent generalized 
quantum Langevin equation (GQLE) in c-numbers2£i2£iSii2£. To achieve this we proceed in two steps. We first carry 
out the quantum mechanical average of Eq.J2J). 



<<?') + f dt'j(t - t'){k{t')) + (V'(q)) = (F(t)) (5) 
Jo 

Where the average (....) is taken over the initial product separable quantum states of the particle and the bath oscil- 
lators at t = 0, |</>){|ai)|a2) I^w)}- Here |</>) denotes any arbitrary initial state of the particle and \di) corresponds 

to the initial coherent state of the i-th bath oscillator. \oti) is given by \cti) — exp(— |ai| 2 /2) ^^° =0 (a™ i /v / "?)|^i) J 
on being expressed in terms of the mean values of the co-ordinate and momentum of the i-th oscillator, (<?i(0)) = 
(\/h/2uji)(ai + a*) and (pi(0)) = i^fkUi/2 (a* — cti), respectively. It is important to note that (F(t)) is a classical-like 
noise term which, in general, is a nonzero number because of the quantum mechanical averaging over the co-ordinate 
and momentum operators of the bath oscillators with respect to initial coherent states and arbitrary initial state of 
the particle and is given by 



(Ht)) = E [{&(°)> - <9(0))K-cos^t+ (^(0))*] 



/2 



siiiLUjt 



(6) 



It is convenient to write Eq.© as follows; 



(q)+ / dt' 7 (t-t , ){q(t')) + {V'(q))=f(t) (7) 



(i 



where we let the quantum mechanical mean value (F(t)) = f(t). We now turn to the ensemble averaging. To realize 
f(t) as an effective c-number noise we now assume that the momentum (pj(0)} and the co-ordinate (qj{0)} — (<?(0)) 
of the bath oscillators are distributed according to a canonical distribution of Gaussian form as, 



so that for any quantum mean value Oj((pj(0)) , {(<?j(0)) — (q(0)}}), the statistical average (....) s is 



{Oj). = / O j P 3 dfe(0)) d{(q 3 (0)) - (q(0))} (9) 



Here fij indicates the average thermal photon number of the j-th oscillator at temperature T and is given by 
Bose-Einstein distribution fij = l/[exp(hujj /kT) — 1] and J\f is the normalization constant. 

The distribution Eq.(JHJ and the definition of the statistical average over quantum mechanical mean values Eq.@ 
imply that f(t) must satisfy 



and 



(/(*)). = (10) 



(f(t)f(t')) s = \ J2 k o fa* ( coth ^f) cos ^(* - f ') 
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or in the continuum limit 



1 f°° f Hlj \ 

(f{t)f(t')) s = - J dLO K (u) pM tkJ ( coth — J cos w(t - 

ee c(t-0 (11) 

That is, c-number noise f(t) is such that it is zero-centered and satisfies standard fluctuation-dissipation relation 
as known in the literature. 

We now add the force term V'((q}) on both sides of the Eq.{7|) and rearrange it to obtain 



q(t) + V'(q) + f dt> 7 (i - 0) = f(t) + Q(t) (12) 
Jo 

where we put (q(t)) — q(t) and (q(t)) — p(t) ; q(t) and p(t) being quantum mechanical mean values and also 



Q{t) = V\q) (V'm (13) 
represents the quantum correction to classical potential. 

Ea. l|12|l offers a simple interpretation. This is that GQLE Ea. (|12|l is governed by a c-number quantum noise f(t) due 
to the heat bath, characterized by the properties Eci. (|10|) . Eq.JTTJ and a quantum correction term Q(t) characteristic 
of nonlinearity of the potential. Our general results contain Q(t) in all orders formally. However for a practical 
calculation we need a recipe for calculation of Q(i). This has been discussed earlier in several contexts2i2£i2£i2L2iLi£i. 
For the present purpose we summarize it as follows: 

Referring to the quantum mechanics of the system in the Heisenberg picture one may write, 



q(t) = (q(t))+Sq(t) 

P(t) = (p(t))+5p(t) (14) 

Sq(t) and 6p(t) are the operators signifying quantum corrections around the corresponding quantum mechanical 
mean values q and p. By construction (Sq) — (Sp) — and [6q, Sp) = ih. Using Eq. i|14|) in (V'(q)) and a Taylor series 
expansion around (q) it is possible to express Q(t) as 



Q(t) = -J2-iV^(q)(6q n (t)) (15) 

n>2 

Here V^ n \q) is the n-th derivative of the potential V(q). To second order Q(t) is given by Q(t) = —\V"{q){5q 2 (t)) 
where q(t) and (Sq 2 (t)) can be obtained as explicit functions of time by solving following set of approximate coupled 
equations Eq.(JTB|) to Eq.(JTHl together with Eq.JHJ) 

(6"% = (Sq5p + dpSq) (16) 
{5qdp + 5pSq)= 2{5p 2 ) -2V"(q){5q 2 ) (17) 
(S :2 p) = -V"(q)(SqSp + 6pSq) (18) 

While the above set of equations provide analytic solutions containing lowest order quantum corrections, the 
successive improvement of Q(t) can be achieved by incorporating higher order contribution due to potential V{q) and 
the higher order effect of dissipation on the quantum correction terms. In Appendix A we have derived the equations 
for quantum corrections upto forth ode»^. Under very special circumstances, it has been possible to include quantum 
effects to all orders22^. The present theory thus takes into account of the anharmonicity as an integral part of the 
treatment. This is somewhat different from usual classical theories where anharmonicity is incorporated from outside 
as a finite barrier corrections^^. We mention here that centroid molecular dynamics simulation methods have been 
developed to calculate the expression for rate constant which is related to flux-flux correlation functions involving two 
nonlinear operators within bilinear system-reservoir coupling^. Apart from anharmonicity of the system potential it 
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is also possible to consider nonlinear coupling, e.g. , ■ 5 Kj /(g) <7j (beyond bilinear) between the system and the 
reservoir. In the present context we, however, envisage two difficulties. First, the potential in which the particle moves 
gets significantly modified because of the absence of any counter term in the Hamiltonian as in Eq. QJ . Secondly, the 
operator Langevin equation of motion would contain multiplicative noise for which a linear dissipation fails to satisfy 
the usual fluctuation-dissipation relation. Because of these reasons we confine ourselves within traditional bilinear 
coupling. The method can, however, be readily extended to anharmonic bath comprising, say, of Morse oscillators or 
of two level systems, for which coherent states are well known (Morse oscillator is amenable to a theoretical description 
within SU(2) Lie algebra and the associated coherent states 51 , while the two-level system in terms of Pauli matrix and 
Radcliffe coherent states^) for constructing quantum Langevin equation since the quantum mechanical averaging with 
these states can be computed exactly in the same way as we proceeded from Eq.J2J to obtain Eq.JSJ. The present 
approach has recently been utilized by us to deriv o 24 ! 26 generalized quantum Kramers' equation for calculation of 
escape rate over the barrier by extending the flux-over-population method to quantum domain. 

B. Conditional probability function of c-numbers at the barrier top 

We now show that the present scheme allows us to derive a quantum conditional probability function in c-numbers 
W(q,p, t; qo,Po, t = 0) for an ensemble of trajectories starting from a specific initial condition q — q , p = p at 
t = and following the quantum Langevin dynamics Ea. (jl2|l . Since we will be primarily concerned with the reactive 
flux at the barrier top located at q — qo{— 0) it is convenient to linearize the potential V{q) at this point as 
V{q) = V(0) — \<~o\q 2 where = V"\ qo —o refers to the oscillator frequency of the inverted well at q = qo. We are 
then lead to the following Langevin equation; 
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q - L0 2 b q + / dt' 7 (i - t') q(t') = f(t) + Q(t) (19) 
Jo 

Laplace transform of Ea. l|19fl yields the formal solutions for q(t) and p{t) as follows; 



q(t) = (q(t)) s + / C v (t-r)f(T)dr (20) 







where 



p(t) - (p(t)) s + f C v (t- r)/(r) dr (21) 
Jo 



(q(t)) s = p Q C v (t) + q C q {t) + G(t) (22) 
(p(t)) , = p C v (t) + q Q C q (t) + G(t) (23) 



G{t)= f C v {t - t) Q(t) dr (24) 
Jo 



ft 

,2 



C,(t) = l+ut C v {t) dr (25) 
Jo 

Here C v (t) is the inverse Laplace transform of C v (p), i.e. 

C v (t)=L- 1 [ M 2 + M 7(/i)-a; b 2 ]" 1 (26) 
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with 



7(A«) = / dt 7 (t) e-"' (27) 



the Laplace transform of 7(t). Several remarks are now in order. 

First, the solutions of Eg . (|20|) and Ea. (|21|l of the Langevin equation Ea. H12() are the quantum mechanical mean 
values distributed around their statistical average (q) s and (p) s . The distribution is due to Gaussian, c- number 
quantum noise f(t) which obeys Ea. (jl()fl and Ea.( |ll|> . 

Second, the averages (q) s and (p) s not only depend on the relaxation functions C v (t) and C q (t) in the same way 
as in the corresponding classical theory^ but also on the convolution integral G(t) which is a quantum contribution 
and takes care of the interplay of dissipation and nonlinearity in the dynamics. It is important to note that although 
we have linearized the potential V(q), Q(t) remains finite and nonzero. This is because linearization involves local 
dynamics around a specific q (say, around q = for the barrier top) which is a quantum mechanical mean value (q) . 
Q(t) as expressed in Ea. (|15|l implies quantum fluctuation around this mean value manifested as shifts in (q) s and 
(p) s as shown in Ea. H22l) and Eci. (|23|) respectively. Making the quantum fluctuations (5q n (t)) or Q(t) equal to zero 
would mean that the quantum mechanical position of the particle is exactly specified. Thus nonzero finite quantum 
fluctuation around q = takes care of uncertainty relation. In other words to consider the conditional probability 
distribution function which originates as a 5-function in position (means quantum mechanical mean) at the barrier 
top, we must have to associate a finite minimum quantum fluctuation or uncertainty around this point. 

Third, since no approximation has been made on Q(t), the solutions Ea. (|20|l and Ea. l|21|l for q{t) andp(t) formally 
take care of quantum effects to all orders. In actual practice, however, they can be calculated order by order to a high 
degree of accuracy. 

Fourth, the quantum effects enter into the dynamics in two different ways. The quantum nature of the system is 
manifested through nonlinearity of the potential while the heat bath imparts quantum noise in a thermal environment. 

f(t) is a c- number Gaussian random process. Since the sum of the random processes is Gaussian, so is the integral 
in Ea. (|2U[l and Eq. l|21[l we conclude that q— (q) s and p — (p) s are Gaussian random processes. Therefore the first and 
the second moments suffice to determine the distribution function. Making use of Eq. I|2(J|) and Eq. I|21|) we construct 

the second moments of q — (q) s ^= J* C v (t — r)/(r) drj and p — (p) s C v (t — r)/(r) dr^j as 

Au.it) = ((q - (q) s ) 2 )s = 2 / C v (h) dh f 1 C v (t 2 ) c(h - t 2 ) dt 2 (28) 



\2\ 



A 22 {t) = ((p - (p) s y)s = 2 / C v {h) dh / C v (t 2 ) c(h - t 2 ) dt 2 (29) 

Jo Jo 



Ai 2 (t)=A 21 (t) = ((q-(q) s )(p-(p) s )) s = 2 C v (h) dh C v (t 2 ) c(h - 1 2 ) dt 2 (30) 

Jo Jo 

We are now in a position to write down the probability distribution function of q, p conditional on q(t = 0) = qo 
and pit = 0) = po 



W(q, P ,t;qo, P o,t = 0) = x exp 



' {A 22 {q-{q) t 



2(detA) 

2A 12 (q-(q) s )ip-(p) s ) + A 11 {p-(p) s ) 2 }} (31) 

Here detA refers to the determinant of the matrix A. The conditional probability W therefore describes the 
transition probability in a c-number quantum phase space. Two remarks are pertinent at this point. First, it is easy 
to check that the W reduces to the classical distribution function in the limit ft — * 0, i.e. , when G(t) — > and 
coth(/iWfc/2fcT) goes over to kT. Second, since W{q,p, t; qo,Po, t — 0) concerns dynamics at the barrier top where 
strong nonequilibrium situation prevails, W does not go over to an equilibrium distribution because of hyperbolic 
nature of the trajectories emanating from q = qo = 0. However, this phase space distribution function serves as a 
basis for calculation of average of the relevant dynamical quantities at the barrier top. 
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C. Equilibrium distribution in the reactant well 

In addition to the distribution at the barrier top W(q,p, i; Qq,Pq, t = 0) it is necessary to consider further the 
equilibrium distribution of the particles in the reactant well. As usual a good description of the kinetics can be 
obtained by linearizing the potential V(q) at the bottom of the well (say at q = ) so that V(q) = V(0) + \u\\q 2 
where ujq corresponds the frequency at the bottom of the reactant well. The Langevin dynamics can be obtained from 
Ea. (|19|l with the replacement of — uj 2 by Wq. The essential changes that are needed for the solution of q(t) and p(t) 



C v {t)=L- 1 C{^) 



(32) 



C g (t) = l-wg / C v {t) dr (33) 
Jo 

The first and second moments of the distribution W can therefore be calculated by employing the relations Eq. i|20|) 
to Eq.||2U and Eq.JSHJ) to Ea.pnfl. Since W(q,p, t; qo,po, t — 0) pertains to the dynamics in the reactant well where the 
trajectories evolve around a stable elliptic fixed point, the distribution W eventually reaches an equilibrium thermal 
distribution in the long time limit, i.e. , 



Ct t ->ooW(q,p,t;q ,po,t = 0) = W eqm (q,p) 



(34) 



Since an equilibrium distribution does not depend on the way in which the final state is reached, one may consider 
the Markovian limit for simplicity but without any loss of generality. In this limit the density distribution of the 
reservoir modes is flat, i.e. , k(lu)p(lu) — 27/71", 7 being dissipation constant. This results in 



and 



7 (f) = 2 7 S(t) 



c(t) = i 7 hwn coth ) 



(35) 



(36) 



from Eq.Q and Ea. l|ll|) . respectively. From Ea. l35fl we have 7(/u) = 7 which when put in Eq. l|32[l gives after inverse 
Laplace transform 



2uj[ 



C v (t) 



where = — ^ + u>[ and /i' 2 = — ^ — with u>'-y = y \ — . By Eq. we obtain 



(37) 



J_( e Mit_ 1 )_i_( e M2t_l) 

Ml A*2 



It is easy to check the stationary limits for the model 



£-tt^<x>C v (i) 



Ctt-+aoC q (t) = 

Ctt^C^t) = 



(38) 



(39) 
(40) 



An important quantity that contains quantum contribution due to the system is G(t) as defined in Ea. H24|) . Putting 
Ea.l|57|l in Eq.lEU) we write 



r Q(t -r)dT- 



Q(t - t) dr 



(41) 
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For the bound system, i.e. , in the neighborhood of the minimum of the well the solutions q(t) and (5q n (t)), in 
general, are oscillating in nature, while the exponentials in Ea. (|41|l are decaying in time. Keeping in mind the nature 



of the system, and expanding Q{t — r) in Taylor series as Q(t — r) = J2 r , 
derivative of Q(t), the long time limit of G(t) is given by 



(-!)'■ 



Q m (t) t" 1 , Q m being the m-th time 



£tt^ooG{t) — - — - 
2uj\ 



(-l) m Q m (t) 



(42) 



Since from Eci. (|15fl the leading order quantum correction is given by Q(t) = —\V"'(q){dq(t) 2 ) it is now easy to 
check that for the time derivative terms in Ea. l|42[l . this correction leads to (and also for higher order correction terms 
in Q(t)), 



CU^Git) = (43) 
for quartic nonlinearity, where we have made use of the relations Eq. I|22l) and Eq. l|23l) . Similarly one may show 



Ctt^&it) = 

By virtue of the relations Eq.JSSJ, Eq.gOJ, En.lj^fy and Ea.ljll}l wc have 

(q) s = o (p) s = o 

Furthermore making use of the relation Eq. 13(j|) for c(t) we obtain from Eq. 
for the second moments. 



(44) 



(45) 



to Eq. (|30|) the following expressions 



An(t) = 7 hujQ coth 



\2kTjJ 



C?(ti) dh 



hujQ \ 



A 22 {t) = 7 ^ coth^^J CUt^dh 

A 12 {t) =A 2 i(t) - ^° cotll (||f) f C v (h) C v (h) dh 
Putting Eq. I|37|l in Eq. (|46|1 to Eq. (|48|l we obtain the long time limits of the second moments as 



Ctt 
Ct t 



A 2 2(t) 

A 12 (t) 



A coth (^M 

2lu \2kT J 
2 C ° \2kT 







(46) 
(47) 
(48) 

(49) 

(50) 
(51) 



Next we use the relations Ea . (|45|l and Ea. (|49|l to Ea. (|51|l in Ea. H31(l to obtain the equilibrium distribution of 
quantum mechanical mean values, q, p as 



W eqm (q,p) = Z 1 exp 



^ + V(q) 



2~ COth ( 2 fcT, 



(52) 



where Z~ x is the normalization constant and V(q) has the linearized form in the neighborhood of the bottom of the 
reactant well characterized by the frequency ljq. The width of the distribution ^p- coth ) reduces to kT in the 
classical limit and one obtains the usual Maxwell-Boltzmann distribution. We thus note that Ea. H52(l has the same 
form as Eq.JHJ) (since, n(x) + 1/2 = \ coth(x/2)). In fact since theq potential at the reactant well is assumed to the 
harmonic, it seems consistent to use the same ansatz that was used for the harmonic bath modes. This has turned 
out to be the case. 

The c-number conditional probability distribution W(q,p, t; qo,Po, t — 0) describing the dynamics at the barrier top 
and the thermal distribution W eqm (q,p) describing equilibrium in the reactant well are the two main results of this 
section. They are primary elements for calculation of statistical averages used in the expression for reactive flux that 
follows in the next section. 
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III. C-NUMBER PHASE SPACE FUNCTION FORMULATION OF REACTIVE FLUX THEORY 



A. A general expression for quantum transmission coefficient 



We now consider q as the reaction co-ordinate and the potential V(q) defines the dividing surface conveniently 
located at the transition state q = qo which separates the reactant state from the product state. While the reactant 
particles are thermally equilibriated deep inside the reactant well strong nonequilibrium situation prevails at q = qo. 
The key element in the reactive flux formulation is that the nonequilibrium concentration relaxes to the equilibrium 
concentration on a timescale at which the equilibrium correlation function decays. Based on a timescale t obeying 
p > t > r c where ^ is the decay time of equilibrium fluctuation and r c is correlation time, the reactive flux expression 
for rate coefficient in the form of Green-Kubo transport coefficient has been derived by various wor kers22i22i2iii&i2£iiiI. 
Our aim in this section is to provide a c-number quantum phase space representation of the rate coefficient. The present 
formulation is a simple quantum translation of the reactive flux formulation advocated by Kohen and Tannor— a few 
years ago. 

In the spirit of classical reactive flux theory the quantum rate coefficient k(t) can be expressed in terms of the 
following correlation function 



k(t) 



(&P (go) dp[q(qo,Po,t)]) qs 

(On(qo))qs 

(PO 5(<l0) &p[q(qO,PO,t)])qs 



<M<zo)) 



(53) 
(54) 



qs 



Here q is the quantum mechanical mean values of reaction co-ordinate and p is the corresponding momentum. 8p [q] 
is a step function implying 1 if q > and otherwise. Noting that R and P refer to the reactant and transition 
states, respectively, we have 6r + 6p = 1. qo, Po correspond to initial co-ordinate (at the barrier top; q = qo = 0) 
and momentum, respectively. The averaging (...) gs is carried over an initial equilibrium distribution, which in our 
case is the c-number quantum distribution W eqm (q,p) of q, p as given by Ea. H52|) . The subscript 'qs' in the average is 
labeled to make a distinction from its classical counterpart. Thus when written explicitly in terms of Ea. (|52|l . Ea. 1541) 
reduces to the following form 



k(t) = 



Z~ 



(0 R (qo)), 



n+oo 

dqo I dpo exp 

-oo 



2 



V(qo) 



^rcoth(^) 



Po 8(q ) P [q(q o ,Po,t)} 



(55) 



For a single realization of a stochastic trajectory the characteristic function 9p[q(qo,po, t)] is either or 1. However 
for an ensemble of trajectories starting from qo, po, 0p[q] must be replaced by an average #p[<?], which may assume any 
value between and 1. In explicit form 6p[q] signifies an integral over the c-number quantum conditional probability 
distribution function W(q, p, t; qo,Po, t — 0) for the ensemble of particles as given by Eg. (|31fl so that 



(■+oo (- + OC 

0p[q(qo,po,t)] = x(P0,t) = / dq dp W(q, p, t; q = 0, p , t = 0) (56) 

JO J -oo 

X is the usual "reactivity index" and takes care of recrossing the transition state. With this replacement of 
9p[q(qo,po, t)] by its average according to Ea. (|56ll we obtain the expression for quantum rate coefficient as 



k(t) = 



(0p(qo))qs 
z- 1 



dqo 



dpo exp 



h^J (no + h) 



po $(qo) x(qo = o,p ,t) 



(0 R (qo)) 



dpopo exp 



qs J - oo 



I + ng = o) 

Tuoo(no + h) 



x(qo = 0,Po,t) 



(57) 



In order to extract out the dynamical contribution to the rate coefficient it is convenient to define the latter as 



k{t) = «(t) k QT 3T 



(58) 
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where n(t) is the quantum transmission coefficient and Icqtst corresponds to the rate coefficient according to 
transition state theory. In the following we outline an expression for kqxsT- 

It has been shown that thermalization of the reactant in the well can be described by the distribution Ea.l)52[l. 
Assuming that this equilibrium distribution holds good even at the barrier top qo Eq. l|52[) reduces to 



W eqm {q =0,p) = Z 1 exp 



\ + V(qo = 0) 



(59) 



and we may calculate the rate coefficient as 



r 

kQTST = I 
JQ 



V W eqm (q = 0,p) dp 



(60) 



which represents the positive average velocity of the particles at the barrier top. The recrossing of the trajectories 
is therefore completely ignored. Explicit integration in Eq.JHOJ) using Eq. l|59|l yields 



kQTST = Z 1 hu (n + -) exp 



V(qo = 0) 
Swo("-o + h) 



(61) 



We now note that the total reactant concentration in the left well is given by equilibrium average over the charac- 
teristic function (9R(qo)) qs which is obtained by 



q +a r + oo 

dq / W eqm (q, p) 

q' — a J — oo 



(62) 



where W eqm (q,p) is given by Ea. H52|) with the linearization of the potential well at q = q' such that V(q) = 
V(q') + ^Wq (q — q') 2 . By setting V(q') = and calculating the normalization constant Z^ 1 = 2Trh(n + |), we have 



7R) qs 



(63) 



Now putting the value of (6ii) qs in Eq.JSTJ) and dividing k(t) by /cqtst as given in Ea. H61|) we obtain the expression 
for quantum time dependent transmission coefficient nit) as; 



K(t) 



1 



huj (n 



-Hoc- 



Po exp 



2 



^wo(«o + I) 



xiqo = 0,p Q ,t) dp 



(64) 



where the fractional reactivity index x is given by Ea. (|56|l . i.e. 



x(po,i)= d q dpW(q,p,t;q = 0,p ,t = Q) 

JO J -oo 

On explicit integration over p using Eq. I|31(l the above expression gives 



(65) 



x(Po,t) 



(2 7 rA U ) 1 /2 



dq exp 



(q-(q)sf 



2A U 



Furthermore with the transformation variables h — g 2y ^| 3 , Eq. (|66|l reduces to 



(66) 



x(Po,t) = —r 



/2 



dh e~ h2 



[2/lll(t)] 1/2 



(67) 
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The expression for n(t) is given by 



K(t) = 



dp po exp 



13. 

2 



Hu} (n + |) 



,1/2 



dh 



(68) 



The primary quantities that determine n(t) are therefore (q(t)) s and An(t) as given by Ea. (|22[) and Eq. 
respectively, pertaining to the dynamics at the barrier top. Putting the value of (q(t)) s in the lower limit of the 
second integral and carrying out integration by parts we obtain the general expression for quantum transmission 
coefficient n(t) as follows; 



«(*) 



C v (t) 



1/2 



exp 



G 2 (t) 



2{C 2 (t) frio (n 



(69) 



The expression Eq. i|69|) is the direct quantum generalization of Kohcn-Tannor formulae^ for classical time dependent 
transmission coefficient. The classical formulae can be recovered in the limit kT fkuo, so that fiwo(^o + 5) an d 
fkOb(nb + \) reduce to kT and G(t) — * 0. A few pertinent points are to be noted: First, since no approximation has 
been made on the order of quantum correction in G(t), the expression takes care of quantum effects to all orders and 
thus retains its validity even in the limit T — > 0, i.e. , in the vacuum regime. Second, since we have not assumed 
any specific form of the density of modes of the bath oscillators so that k{u)p{u) is arbitrary, the theory developed 
so far is equipped to deal with arbitrary noise correlation. Third, it is important to note that important quantities 
in the expression Ea. l69fl . C v (t), Au(t), and G(t) are the key elements required for calculation of the dynamics 
of transmission coefficient n(t). While C v (t) is the same as in the classical theory, Au(t) and G(t) are quantum 
mechanical objects. An(t) carries the imprints of quantum fluctuation-dissipation relationship characterizing the 
quantum nature of the heat bath. G(t) on the other hand, as mentioned earlier is a reflection of the characteristic 
of the nonlinear potential. The interplay of the dissipation with the nonlinearity of the system is expressed through 
this convolution integral. A complete calculation of the transmission coefficient thus requires the explicit knowledge 
of these three quantities. 

In order to proceed further it is necessary to specify the nature of the time dependent friction j(t). This is crucially 
dependent on the distribution of the density of modes of the bath oscillators. In the present case we assume a specific 
form of Lorentian type, i.e. ; 



k(ui)p(u>) 



A 



7T V 1 + LO 2 T? 



(70) 



where A is the damping constant in the Markovian limit and r c is the correlation time. This form of density by 
virtue of Eq.QJ results in, 



7(«) 



.4 



-t/T c 



A a e 



(71) 



where a = l/r c . 

Depending on r c we consider two different limits; Markovian and non-Markovian friction. We consider the former 
case first. 



B. Markovian friction 



For a — > 00 or r c — > 0, we have k(w)p(lu) = 2A/tt] j(t) and its Laplace transform thus reduce to 



7(i) = 2 A 5(t) 
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and 



7(/x) = A 

respectively. According to Ea. l(2*Hjl C v {t) is then given by 



(72) 



C v (t) 



2wi 



(73) 



where /ii^ = —4 ± and lox 



A 1 



Furthermore the correlation function as given by Eq. I|ll|) yields 



c(ti - £2) = ^ ^4 coth (^f) <*(*i - 
Making use of the last two relations for C v (t) and c(ti — t^) for the barrier top in Eq. l|28|l we obtain 



(74) 



Ml{t) = ^ huj b (n b + ~) 



1 

2^1 



2/^2 



.4 



(75) 



Having known C v (t) and j4n(£) in the Markovian limit we are in a position to use them in the general expression 
for quantum transmission coefficient Ea. (|69|l to obtain, 



K(t) 



x exp 



G 2 (i) 



_J^ e 2/ii« j ]_ e 2fi 2 t i 2_ e -Ai _ 2 ^i 



1/2 



^o(n + ^) 



2w? 



--4*) _|_ 2AhuJ b(n b + ^) 



_l_ e 2/iii i _J_g2/j 2 t + 2. e ^ At — 2 " 1 
2/ii 2^2 ^4 



(76) 



It is instructive to work out the assymptotic limit of K,(t). For a bistable potential of the form V(q) = aq — bq 2 , the 
leading order quantum correction is given by Q(t) = —12 a q(t) (Sq 2 (t)). Since the dynamics pertains to the barrier 
top we have approximately [see Appendix B] 



G(t) = - 



3« 



Ldb Wl A 2 

which may be put in the form 



h 3 u (n + -) 



1/2 



— sinh2w b i (e" 1 * + e^ 2 *) 



2 cosh 2 u b t (e^-e" 21 ) 



(77) 



G(<) = 



3 a 



A 2 



H 3 u; (n Q + -) 



1/2 



2 2 



in the long time for which lnwi <C 2uj b t and lnwb <C 2wbi. Eq.J78J) then gives 



(78) 



G(t) 



3 a 



A 2 



H 3 uJo(no + -) 



-,1/2 



(79) 



Putting Eq. l|79|l in Eq. l|76|l and using the limit t — > oo we obtain the assymptotic expression for the transmission 
coefficient in the Markovian limit. 



k(oo) = 



^ A huj b (n b + ^) 
^1 Rw (n + i) 



exp 



18 a 2 h 2 fix uj Q (n Q + i) 



A 4 lo 2 (/ii u (™o + 5) + ^ ^b(wfc + 5)) 



(80) 
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In the classical limit as h — » oo, hu>b(jib + h) and Hu)o(no + 5) both goes over to kT; Eq. l|80|l results in 



k(oo) 



1 + ^ 

Ml 



1/2 



(81) 



Since /ii + /12 = — ^ and /ii/12 = — w 2 the above relation can be put into the following form 



, v Ml 

K(OO) = 



(82) 



The classical Markovian coefficient, i.e. , the ratio of the positive unstable root /ii to the frequency of the barrier 
LOb is therefore reduced to Kramers' prefactor. 



C. Non-Markovian friction 



We now return to the Eq.(7U|) and Eq.{7TJ- Corresponding to Eq.(7TJ we have = from Eq.(|2Z|). C v (t) is 
therefore is given by 



C v {t) = L- 1 



2 2 . V Aa 
M ~ u b + — ■ — 

M + a 



(83) 



Following Kohen and Tannor— we consider the case when a is small, i.e. , the relevant zeros of the quantity in the 
parenthesis in Ea. H83(l may be written as the solutions of 



- uj 2 + A a = 



(84) 



Depending on the relative magnitudes of uj 2 and Aa it is convenient to discuss two situations ui 2 — Aa > and 
ui 2 — A a < separately. 
(i)wj-ia>0 

Here the two roots become /11.2 = ±y/uj 2 ~ A a, with /«i = —[12 and C v (t) is given by 

sinh jixt 



C v (t) = 



Mi 



Since k(ui)p(u)) — we have from Eg. Hll|) the correlation function as follows 



. 1 Z" 00 , /2A\ a 2 ^ 
c(ti —h) = -zl da; — — - ftwcoth 



V 2kT 



COSU>(ti — t-z) 



Putting Eq.JSBJ and Eq.JSSJ in Eq. <|2*5|) it is possible to express Au(t) as 



A a 2 H sinh 2 ui t , 
A n (t) = , — — M - 



2 A a 2 ftsinhuit 2 ^4 a 2 ficoshuii ^ A a 2 H 
— N — CM 

TTfll 7T 7T 



where 



M = 

N = 

O = 

P = 



id 



/•OC 



(a 2 


+ uj 2 ) (to 2 + 


Mi) 


cj 2 


sin ut coth 


2kT 


(a 2 


+ U) 2 ) (U! 2 + 


Mi) 2 


w 


cos cot coth 


2kT 


(a 2 


+ UJ 2 ) {uj 2 + 

uj coth^ 


Mi) 2 



/„ (a 2 +w 2 ) (w 2 + /i 2 ) 5 



(85) 

(86) 

(87) 

(88) 
(89) 
(90) 
(91) 
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The nonlinear term G(t) in Eq. (|69|l can be obtained as before by considering the leading order quantum correction 
as shown in Appendix B for the case — Aa > 0. This is given by 



G(t) 



3 a 



2 u>b Mi A a 



, ■ , 1 

n gjq [no + - 



iV2 r 



Hi smh(/zii) H cosh(/iii) smh(— — ) 

uji, A 



0Jb smh( — — J cosh(/iit — J — /ii cosh(— — ) smn(/iif —J 

yd -Tl J\ J\ 



2m 



,2cOb, 



2/xi 



(92) 



Having obtained C v (t), An(t) and G(i) from Ea. (|85|l . Ea. (|87|l and Ea. 192(1 . respectively, we may write down the 
quantum transmission coefficient using Ea. l|(j9fl in the adiabatic regime (w 2 > Aa) 



K(t) 



rsinh_££i£i 

i Hx -I 



lh 2 n\t 



Aa 2 h 



exp 



p£ N-2cosh/j, 1 t + 2P 



1/2 



{^ (n + l) 2 ^ 



; x A^B^ smhVit M _ 2 B inh Ml t N _ 2 COsh Ml i O + 2P)} 



(93) 



In order to check the workability of the above method we first test it against two typical benchmark calculations. 
To this end we first employ the double-well potential as used earlier in exact numerical calculation of quantum 
rate coefficient in Ref.[22] with scaled parameters (actual parameters in Table I. of DW1 in Ref.[22]) Ej, = 26.188, 
Lu b = 1.0, uj = 1.414 for two different scaled temperatures kT = 2.617(300^) and kT = 1.744(200isT). The results for 
assymtotic transmission coefficient as a function of damping constant from Fig. (9) of Ref.[22] are compared against 
our analytical expression for the same according to Ea. (|93|l in Fig.(l). The agreement is excellent at relatively higher 
temperature i.e. at kT — 2.617(300if). However it is apparent that at lower temperature the analytical results is 
slightly higher than than exact one. Keeping in view of the fact that the maximum deviation is less than 8 percent 
and that the analytical quantum correction in Eq. I|93|I is of lowest order, the result is fairly satisfactory. For a 
more accurate calculation one would require higher order corrections in Q(t) which can be obtained from numerical 
solution of the equations in Appendix A and be included in G(t) to be used in Eq. l|63|l . In what follows we also 
make use of the classical results of Cohen and Tanno»2£ for a consistency check of our quantum calculation in the 
approximate limits as we proceed. We now show the detailed behavior of n(t) as a function of time for different values 
of temperature T, damping constant A and correlation time r c (= 1/a). To this end we first compare the quantum 
nature of transmission coefficient as a function of time against its classical counterpart for A — 50, t c = 100, kT = 1.0, 
tub = f -0, ojq — \[2 and a — 0.01 in Fig. (2). It is apparent that there is a significant increase in quantum values. 
The variation of temperature has been depicted further in Fig. (3). At higher temperature the quantum coefficient 
merges into the classical coefficient of Kohen and Tannor 36 , An important new content of the present theory is 
that it is equipped to deal with zero point fluctuations. Around this teperature the transmission coefficient hardly 
departs from its maximum value. In Fig. (4) we show the explicit variation of assymptotic transmission coefficient 
with temperature. While classically n(t) is independent of temperature, quantum nature of the heat bath imparts a 
temperature dependence which is particularly strong in the low temperature regime. This is one of the most distinctive 
features of the present calculation. 

Fig. (5) shows the two sets of transmission coefficient n(t) as a function of time for different values of dissipation 
parameters A ranging from 30 to 90 one for quantum and other for the classical. The input parameters are kT = 1.0, 
u>b = 1.0, ujq = v2, a = 0.01 and r c = 100. It is apparent that there is marginal increase of the corresponding quantum 
values when compared to the classical transmission coefficient. In Fig. (6) the assymptotic transmission coefficient is 
plotted as a function of damping constant. The typical behavior in the intermediate to strong dissipation regime is 
observed which is in complete conformity with the earlier results. Fig. (7) shows an illustrative variation of the time 
dependent quantum transmission coefficient for several values of correlation time r c . With increase in t c n{t) increases 
significantly - a trend which is similar to what is observed in the classical case. With increase in time the transmission 
coefficient n(t) reaches its plateau value in the Markovian as well as non-adiabatic non-Markovian regime. 

(ii) ul- Aa <0 

Here from the equation for /i, i.e. , /i 2 + {A a ~ uj%) — we obtain two purely imaginary roots //2,3 



fj, 2 ,3 = ± «M2 
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where fa = \J Act — uj 2 . C v (t) in this case is given by 



C v {t) 



sin fat 

fa 



(94) 



Proceeding as before we arrive at the quantum transmission coefficient in the so called caging regime as follows 



«(*) 



r sin ftgt i 

fJ-2 ' 



Aa 2 H 



sin 2 ubf 

fi2 1 Trhu>o{no + ^) 



sj^pt Mi _ jv x _ 2 COS fat Ox + 2Pi 



1/2 



x exp 



G{tf 



2{fi^ (no + i)^f^ 
Here G{t) is given by [see Appendix B] 



7T V ^ 1 



2 sin 
£2 



iVi -2cosfat Ox + 2Pi)| 



(95) 



G (t) = - 3 ° - 

2W6 



H 3 ujq n + 



— fl2 cosh I — — 



sm I /i 2 i - — 



Wf, smh I — — I cos I ^~ 

- fa sin (£t2i) cos (fat) sinh ( — 

LU b \ A 



(96) 



where 



Mi 
Nx 
Ox 
Pi 

The behavior of the transmission coefficient is very much similar to that for the classical case as discussed in detail 
by Kohen and Tannor— . In order to allow ourselves a fair comparison with this theory we calculate according to 
Ea. (|9"5|l and the corresponding classical case for the parameters A — 150, r c = 100, kT = 1.0, LJb = 1.0, uJq = v2 and 
a = 0.01. The results are shown in Fig. (8). It is apparent that the oscillating nature persists even in the quantum 
case although this is accompanied by a loss of amplitude compared to the classical case. In Fig. (9) and Fig. (10) 
we illustrate the dependence of transmission coefficient with temperature and damping, respectively. Apart from 
temperature dependence no significant difference from classical behaviour is observed. 



(a 2 
u? 


+ U) 2 ) (uj 2 - 
sxriLot coth 


Ml) 

huj 

2kT 


(a 2 
u> 


+ w 2 ) {u> 2 - 
coswi coth 


ti? 

2kT 


(a 2 


+ U! 2 ) (U) 2 - 

uj coth^ 


ti) 2 



(a 2 +uj 2 ) (uj 2 -fi 2 2 ) 2 



(97) 
(98) 
(99) 
(100) 



IV. CONCLUSION 



In this paper we have derived a general time-dependent quantum Kramers'-Grote-Hynes rate constant for a chemical 
reaction in a condensed medium. This is based on quantum phase space probability distribution functions proposed 
recently by us to describe quantum Brownian motion in terms of a generalized quantum Langevin equation. We have 
correlated the time dependence of the transmission coefficient in detail with the variation of dissipation, temperature 
and correlation time of the noise. Our results can be summarized as follows: 

(i) The quantum KGH transmission coefficient Ea. H69|l . the central result of this paper reduces to to its well known 
classical limit when h — ► 0. The expression retains its full validity even in the vacuum limit, i.e. , T — > 0. 

(ii) We have worked out in detail the specific form of this transmission coefficient in the Markovian as well as 
in the non-Markovian case. Depending on the dissipation, correlation time and frequency at the barrier top, both 
adiabatic and non-adiabatic situations have been considered with a comparison with the classical theory of Kohen 
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and Tannor 36 . While classical transmission coefficient is independent of temperature its quantum counterpart depends 
on temperature significantly, particularly in the low temperature regime marked by quantum effects. 

(iii) The general expression for quantum KGH transmission coefficient contains quantum corrections to all orders. 
We have shown that the systematic corrections can be obtained explicitly order by order by simply solving a number 
of coupled ordinary differential equations. 

(iv) The expression for the quantum KGH transmission coefficient is a natural extension of classical theory to 
quantum domain and provides a unified description of thermally activated process and tunneling. 

Before conclusion we now make a few remarks on the present method. 

First, since the advent of quantum mechanics several methods of quantization have been in use for description of 
quantum Brownian motion. Of these, Lagrangian and Hamiltonian formulation have been proposed in the frame- 
work of Feynman (path integral) quantizationi^S^ and been most widely used since early eighties. Dynamical 
quantizationi§*2Li& have been tried with some preliminary success. The present method is based on canonical 
Heisenberg quantization procedure and has been recently utilized to propose non-Markovian quantum version of 
Kramers242&, Fokker- Planck^, SmoluchowskiSl and diffusion equations^. The main advantage of this approach is 
that it allows the classical method of solution of equation for probability distribution functions to be applied in the 
same way as one makes use of fiux-over-population method or reactive flux method for calculation of rate coefficients. 

Second, since the phase space function formulation of classical reactive flux theory by Kohen and Tannor is partic- 
ularly useful for visualization of the behaviour of phase space trajectories, a direct extension of this method to the 
quantum domain as done in the present formulation allows us to achieve a similar objective in terms of c-number 
trajectories. It is also important to point out that path integral approach also gives rise to several approximations 
that have a classical flavour^ 3 ^. 

Third, an important advantage of this approach is that one not only obtains the asymptotic quantum rate coefficient 
but also the details of its behaviour at all time analytically including its approach towards the steady state. In this 
sense it is complimentary to the exact numerical work of Ref . [22] . 

Fourth, although the present method employs system-reservoir oscillator model for evolution of quantum mechanical 
rate, in the phase space distribution function approach bath co-ordinates do not appear explicitly in the ultimate 
calculation. One can also avoid multidimensional integration employed in more powerful path integral approach. 
The present method instead takes care of improvement of quantum corrections successively order by order by solving 
coupled ordinary differential equations of Appendix (A). 

In the view of the above discussions we anticipate that the present approach can readily be adopted to numerical 
scheme, for example, to that of classical method of Berne and co-workers^S^S, for calculation of time- dependent quan- 
tum transmission coefficient beyond linearization over the entire range of friction, noise correlation and temperature, 
to explain quantum turnover and other interesting features. The work in this direction is in progress. 
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APPENDIX A: EVOLUTION EQUATIONS FOR HIGHER-ORDER QUANTUM CORRECTIONS FOR 

ANHARMONIC POTENTIAL 

The equations upto fourth order for quantum corrections (corresponding to the contribution of anharmonicity of the 
potential) are listed below. 
Equations for the second order are: 



Acknowledgments 




(SqSp + SpSq), 



2T(Sp 2 ) - V"(SqSp + SpSq) - V'"(5qSpSq), 



(Al) 



— (SqSp + SpSq) 



T(SqSp + 5pSq)2(Sp 2 ) - 2V"(5f) - V"'(6q 3 ), 



Those for the third order are: 
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j t (6f) = 3(5qSpSq), 

j t {5f) = -3T(Sp 3 ) - 3V"(SpSqSp) + V" (~{6?){8?) - \(5p8q 2 5p) + h 2 

-(SqSpSq) = -T(SqSpSq) + 2(S P SqSp) - V" (5q 3 ) - — [(Sq 4 ) - (Sq 2 )) , (A2) 

d (SpSqSp) = -2T(SpSqSp) + (Sp 3 ) ~ 2V"(SqSp5q) 



dt 



ylll 

— ((5q 2 }(SqS P + SpSq) - (Sq 3 Sp + 5 P Sq 3 )) , 



And lastly, the fourth order equations are 

d 



^-(Sq 4 ) = 2(6q 3 5p + Sp6q 3 }, 
dt 

4(^ 4 ) = -4T(Sp 4 ) - 2V"(SqSf + SfSq) + 2V"'{8q 2 ){5p 3 ) 1 
dt 

^-(SfSp + SpSf) = -T(5q 3 Sp + 5p5q 3 ) - 2V"(6q 4 ) - 3h 2 + 6(5p5q 2 Sp) 
dt 

+ V"'(5q 2 )(5q 3 ), (A3) 

^(SqSp 3 + Sp 3 Sq) = ~3T(SqSp 3 + Sp 3 Sq) +2(5p 4 } +3V"(h 2 -2(SpSq 2 5p)) 
dt 

+ W"{5q 2 ){5 P 5q5p), 

— (SpSfSp) = -2T{5 P Sq 2 6p) - V"(Sq 3 Sp + 5p5q 3 ) + {Sp 3 5q + SqSp 3 ) 

+ V'"(6q 2 )(6q5p5q). 



dt 



The derivatives of V(q), i.e., V" or V'" etc. in the above expressions are functions of q the dynamics of which is given 
by Eq. lQ . 

APPENDIX B: CALCULATION OF G(t) 

It is already been pointed out that G(t) as defined in Eg. 1241) describes an interplay of dissipation and nonlinearity 
in terms of the following convolution integral 

G{t) = f C v {t~r) Q(r)dT (Bl) 
Jo 

where Q(t) is given by Ea. p5|l and the relaxation function C v (t) refers to Ea. (|26[l pertaining to the dynamics at the 
barrier top. When Q(t) is taken in full G(t) contains quantum corrections to all orders and is formally exact. However 
for practical calculation we evaluate this quantity order by order. Leading order quantum correction in Eq. l|15|l as 
given by Q(t) — —\V"'(q)(8q(t) 2 ) for the bistable potential V(q) — a q 4 — b q 2 results in 

Q{t) = -12 a q(t) (Sq{t) 2 ) (B2) 

Explicit solution of (Sq(t) 2 ) at the barrier top 24,39 at q — qo = is given by (Sq(t) 2 ) = (5q 2 )t=o cosh2oj6t + 
(Sq Sp+Spjq) t =o s i n ] 1 2a;6t. Furthermore q(t) is known from Ea. H2U|) upto a leading order so that we have 

q(t)= Po C v {t) (B3) 
Since with minimum uncertainty product stat o 24 ' 3 ? we require 

(Sq Sp + SpSq) t=0 = (B4) 
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and so we have 



(Sq(t) 2 ) = — cosh2uj b t 



2ui b 



(B5) 



(B6) 



For an initially thermally distributed po we may write 



Po 



-l 1/2 



huj (n + -) 



(B7) 



With Eq.jBZjl, Eq.jBBj) and Eq.(|B3), Q(t) yields 



Q(t) 



6 a h 3 uj (n + i) 

LU b 



C v (t) cosh.2ui b t 



(B8) 



Putting Eq.jBBJl in Eq. (|Bll) we obtain the general expression for G(t) upto the leading order in quantum correction 



as 



G(t) = -^ 

LU b 



h 3 uj (n + ~) 



1/2 



C w (i — r) C„ (r) cosh 2u b T dr 



(B9) 



The rest of the treatment depends on specific nature of C v (t). 

Case(I) Markovian friction: For this case C v (t) is given by Eq. l|73|l . Explicit evaluation of the integral in Ea. (|B9 
results 



G(t) 



3 a 



OJl UJb 



A 2 



h 3 Lo (no + -) 



1/2 



— sinh2w 6 i (e^ 1 + e" 2t ) - 2 cosh 2 u b t (e^ 1 * - (e^ 2 *) 

UJ b 



(BIO) 



Case(II) non-Markovian friction: 

Similarly the evaluation of the integral in Eq. (|B9|) using C v (t) from Eq. I|85|l for the non-adiabatic case u>1 — Aa > 
leads to 



G(t) = 



3 a 
2 ui b n\ A a 



/ 1MV2 
h 3 ' ' 



3 w ( "o + 2 



/Lti smh(/zii) H 

— smn( — — ) cosh(/iir —J — 



cosh(/iit) sinh( — -) 



//i cosh(— — j smh(//i£ —J 



(Bll) 



where to avoid divergence we have used a cut off for the integral at time t = 1/A. This is a good approximation in 
view of the fact that quantum correction in Eq. (|Bll|l has been considered to a leading order. 



For adiabatic or so called caging regime (w 2 — Aa < 0), on the other hand C v (t) is given by Ea. (|94|) . The 
corresponding G(t) can be calculated to obtain from Ea. (|B9|) as 



G(t) 



3 a 



2 uj b \x\ A a 



h 3 u>o I no + 



1 



-1 1/2 r 



sinh ( — — ^ j cosh {[Lit -^ 2 - 



j) 2 co-,1) ( ^^- \ sinh Ufat — + M2 sinh(// 2 <) - -^cosh(/2 2 <) sinh {^~^\ 



(B12) 
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FIG. 1: Asymptotic quantum transmission coefficient n s t is plotted as a function of damping strength (A) at two different 
temperatures (a) kT = 2.617(30(1R') and (b) kT — 1.744(200if) for the parameter set mentioned in the text. The numerical 
results of Fig. 9 of Ref.(22) are shown by solid circles and squares and compared to our analytical results (solid triangles) 
corresponding to Eq.(93). 
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FIG. 2: The quantum (continuous line) transmission coefficient is plotted against time and compared to classical (dotted line) 
one for the parameter values A — 50, r c = 100, kT — 1.0, oji = 1.0, u>o = V2, a = 0.01. 
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FIG. 3: The quantum transmission coefficient is plotted against time for several values of temperature kT — 0.0 (continuous 
line); kT = 0.3 (dashed line); kT = 25.0 (dotted line) for A = 50, r c = 100, kT = 1.0, co b = 1.0, lo = y/2, a = 0.01. 
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FIG. 4: The assymptotic quantum transmission coefficient n s t is plotted against temperature kT for the set of parameters as 
in Fig.2. 
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FIG. 5: The quantum (continuous lines) transmission coefficient n(t) is plotted against time for different values of damping 
constant (a) A — 30, (b) A — 70, (c) A = 90 and compared with classical (dotted lines) ones for the parameters r c = 100, 
kT = 1.0, co b = 1.0, ujo = V2,a = 0.01. 
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FIG. 6: The assymptotic quantum transmission coefficient n(t) is plotted against damping constant A for the parameter values 
as in Fig.4. 
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FIG. 7: The quantum transmission coefficient is plotted against time for different values of correlation time, r c = 65 (continuous 
line), t c — 75 (dotted line), r c = 95 (dashed and dotted line), r c = 110 (dashed line) for other parameters as in Fig. 4. 
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FIG. 8: The quantum (continuous line) transmission coefficient is plotted against time and compared with classical (dotted 
line) one for the set of parameters A = 150, r c = 100, kT = 1.0, LOb = 1-0, ujo = V%, a = 0.01 (in the caging regime). 



